function Policy=Policy_Approximate(PP,SS,GridPolicy,State)

N_State     =   size(State,1);
CoefMat     =   SS.PolApp.V.BasMat_0\reshape(GridPolicy,[SS.PolApp.V.MatSize_State]);
Policy      =   zeros(N_State,1);

for ii_a=1:PP.ExoState.TempTFP.N
    TempInd     =   (State(:,end)==ii_a);
    if any(TempInd)
        Policy(TempInd) ...
                    =   funbas(SS.PolApp.V.Space,State(TempInd,1:end-1)) ...
                        *CoefMat(:,ii_a);
    end
end